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ABSTRACT 

We measure the fraction of luminous galaxies in pairs at projected separa- 
tions of 5 — 20kpc out to z = 1.2 in the COSMOS field using ACS images and 
photometric redshifts derived from an extensive multiwavelength dataset. Anal- 
ysis of a complete sample of 106,188 galaxies more luminous than My = —19.8 
(~ Ly) in the redshift range 0.1 < z < 1.2 yields 1,749 galaxy pairs. These 
data are supplemented by a local {z = — 0.1) value for the galaxy pair fraction 
derived from the Sloan Digital Sky Survey (SDSS). After statistically correcting 
the COSMOS pair sample for chance line-of-sight superpositions, the evolution 
in the pair fraction is fit by a power law oc (1 + If this strongly 

evolving pair fraction continues out to higher redshift, ~ 50% of all luminous 
galaxies at z ~ 2 are in close pairs. This clearly signifies that galaxy mergers are 
a very significant and possibly dominant mechanism for galaxy evolution during 
the epoch of galaxy formation at 2; = 1 to 3. 
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Subject headings: cosmology: observations — cosmology: large scale strutcure of 
universe — galaxies: formation — galaxies: evolution — galaxies: interactions 
— surveys 



INTRODUCTION 



Galaxy mergers play an important role in galaxy formation and evolution. Not only do 
they change the number density of galaxies in the universe and provide a method for mass 
accretion, but they also significantly alter the shapes and morphologies of galaxies, induce 
bursts of star formation, and provide material for accretion onto black holes. Current models 
of galaxy formation sugg est that galaxy mergers must have been much more common in the 
past than they are now (IToomrd 119771 : iBarnes fc HernquistI 119921 ). 



Though the significant role of merging has long been recognized, quantitative estimates 
of merging rates have varied greatly. Derived merger rates depend strongly on the redshift 
range surveyed, the wavelength band used for imaging, and the limiting magnitude of the 
survey. In addition, all previous studies at z > 0.1 have involved relatively small galaxy sam- 



use pair fraction as a surrogate for the merger rate (e.g., Zepf & Koo 


1989: Burkev et al. 


I994I: Woods et al. 


I995I: Yee & I 


]Uineson 19951: 


Neusc 


laefer et al. 1997 


: Patton et al. 2002: 


Wu & Keel 19981: ( 


ZJarlbere et al. 


2OO0I: Le Fevre et al. 


2OO0I: Lin et al. 


2004): however, the 


pair fraction may not translate so simply (e.g.. 


Abraham 1999). Instead, it samples ealaxv 



pairs at different stages of interaction and assumes that that they will all eventually merge. 
Although most close galaxy pairs may ultimately merge, derivation of the actual merger rate 
from the pair fraction requires an understanding of the merging timescale which is likely to 
be strongly dependent on both the observed separation and the assumed orbital eccentricity. 

Several recent studies use morphological criteria to select disturbed systems in the ad- 



vanced stages of merging (e.g.. lLe Fevre et al.ll2000l : IConselice et al.ll2003l : iLavery et al.ll2004 



Lotz et al.ll2006l ). While this method is more likely to yield objects which have just recently 
merged or are very likely to merge, it has its own deficiencies: information on the progen- 
itor galaxies is often lost and the morphological selection can be very redshift-dependent 
due to surface brightness dimming and morphological K-corrections. The morphological 
appearance of galaxies can be quite dependent on the observed wavelength, particularly as 
one observes more and more of the rest frame UV wavelengths. The galaxy morphologies 
then become increasingly dominated by irregular areas of recent star formation, rather than 
galactic mass structures. When the star formation is particularly intense and clumped, the 
blue-UV morphology may then mimic the disturbed appearance of interacting systems. 
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Given the differences and problems listed above, merhaps it is not surprising then that 
both techniques - pair fraction and morphology - have yielded a very broad range of power- 
laws ((1 + z)"^) for the merger rate - from m = (no evolution) to m ~ 5 (strong evolution) 
as reported in the literature (see §5.2). 



The Cosmic Evolution Survey (COSMOS: IScoville et al.ll2007bl ) offers a new and unique 



opportunity to study the galaxy merger rate versus redshift. The availability of HST imaging 
and galaxy catalogs for the full 2 square degree COSMOS field (Scoville et al. 2007c, Koeke- 
moer et al. 2007, Leauthaud et al. 2007) provides an enormous sample of galaxies, while 
extensive ground-based follow-up imagi ng of the full field in t he optical and near-infrared 



provides accurate photometric redshifts (IMobasher et al.l 120071 ). Here we make use of these 



data to investigate the evolution of the galaxy pair fraction, taking advantage of the high 
sensitivity and redshift accuracy of the COSMOS catalogs. Specifically we use a sample of 
106,188 > L* galaxies to derive the fraction of close galaxy pairs (at 5 — 20kpc separation) 
out to z = 1.2. 

Most earlier studies have used spectroscopic redshifts to confirm that galaxies close in 
position on the sky are indeed at similar redshifts and therefore possibly gravitationally 
bound. These studies have typically probed projected separations less than a few tens 
of kiloparsecs with AV < 1000 km/s. However, the requirement for spectroscopic redshift 
confirmation usually results in very limited sample sizes, and since the full 3-d velocity is still 
not known even for the spectroscopic samples, one still doesn't know that the pair is indeed 
bound. Thus the major utility of the spectroscopic redshifts is to greatly reduce line-of-sight 
contamination or chance superposition. However, since the line of sight contamination can 
be statistically estimated, corrected for, and kept at a low level (by avoiding very large 
impact parameters), we make use of photometric redshifts; these have lower accuracy than 
spectroscopic redshifts but yield vastly larger samples. The photometric redshift sample has 
the added benefit that it is also free of the fiber collision problems and incompleteness which 
often plague spectroscopic samples. 

The structure of this paper is as follows. The data sets used for this project are described 
in §2 and the selection criteria used to detect the pairs are discussed in §3. Our results are 
presented in §4 and discussed in §5 where we compare results with those of previous studies 
on the pair fraction and merger rate evolution. Our conclusions are presented in §6. The 
cosmological parameters used throughout this paper are h = 0.7, = 0.7, and Qm = 0.3. 
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THE DATA 



The imaKiiiK data used for this study were obtained as part of the HST-COSMOS 
project (IScoville et al.ll2007d ). COSMOS originated as an ifS" T Treasury program imaging 
an ~2 0° equatorial field with the Advanced Camera for Surveys (ACS), using the F814W 
filter (I-band). This is the largest contiguous field ever observed by HST. The field size 
was chosen to sample the broad range of large scale structure (LSS) and thereby re duce the 
problem of cosmic variance that plagues smaller survey fields ( IScoville et al.ll2007bl ). Details 
of the ACS images, including their calibration and reduction, are given in Koekemoer et al. 
(2007). In addition to the ACS images, we used the ACS source catalog that is complete to 
/ < 26 Jleauthaud et al.ll2007l l 



The COSMOS project has also obtained extensive multiwavelength ground and space- 
based observatio ns covering the entire f ield. A complete description of the COSMOS data sets 
can be found in IScoville et al.l (l2007bl ). A ground-based photometry catalog was produced 
using imaging data taken with Subaru-SuprimeCam [BVg'Ri'z' NB816), CFHT-Megacam 
(M*,r), CTIO-ISPI (Ks band), NOAO Kitt Peak-4m {K^ band) and SDSS {u' g'r'i' z'). De- 

tails of the observa tions , data reduction and c atalogs can be found in Aussel et al. (2007), 

Capak et al.l (120071 ) , and lTaniguchi et al.l (120071 ) . Using the ground-based photometry, iMobasher et al 



(120071 ) calculated photometric redshifts for approximately one million galaxies within the 
COSMOS field. These photometric redshifts have an accuracy of dz/{l + z) ^ 0.03 out to 
a redshift oi z = 1.2 (IMobasher et al.ll2007l ). The photometric redshift catalog is crucial 
for our study, enabling a sorting of the galaxies in redshift bins to reduce contamination by 
chance, line-of-sight pairs and hence to yield a reliable sample of galaxy pairs. 

At low redshift {z < 0.2), the COSMOS field encompasses too small a volume to provide 
a meaningful esitimate of the local galaxy pair fraction. Therefore, to estimate the "local" 
(2; = — 0.1) fraction of galaxy pair s, we make use of th e published catalog of Sloan Digital 
Sky Survey (SDSS) galaxy pairs by lAUam et al.l (120041 ). This catalog samples a volume in 
the local universe nearly as large as that in the COSMOS field aX z = 1. 



3. ANALYSIS 



3.1. Galaxy Sample Selection Criteria 



The primary criterion for galaxy selection w as an absolute lumin osity cutoff correspond- 



ing to a galaxy, specifically My = —19.8 (IBlanton et al.ll200ll ). This was done partly 



for practical reasons in that galaxies above this luminosity have apparent I band magnitudes 
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of < 26 and are still easily visible out to the redshift limit of this study, z = 1.2, thus 
eliminating the need to correct for incompleteness at higher redshifts. More importantly, it 
is likely that pairs of galaxies brighter than Ly will become components of approximately 
equal mass mergers since most of these galaxies are near the "knee" of the luminosity func- 
tion (since above the "knee", the number density of galaxies drops steeply). Lastly, we note 
that mergers between galaxies with luminosities > L* are likely to be responsible for some 
of the most luminous an d rapid phases in galaxy evolution, includi ng ultraluminous infrared 
galaxies (UL IRGs: e.g., Kim fc SandersI IQQsl: Dasyra et al. 2006 ) and quasi-stellar objects 



(QSOs: e.g.. lVeilleux et al.lbood : iGuyon et al.ll2006l ). To ensure that all the galaxies in our 



sample have reliable photometric redshifts we required that each be detected in at least four 
bands, including band to guarantee a wide wavelength coverage. Applying the above se- 
lection criteria to the 438,356 entries in the ground-based catalog yielded a sample of 106,188 
galaxies > Ly across the entire COSMOS field. Selecting all galaxies with z < 1.2 yields a 
final sample of 59,221 galaxies used to search for pairs. 

For studies which examine the pair fraction as a function of redshift it is critical to 
impose a lower limit to the galaxies' luminosities (or even better their masses) rather than 
an apparent magnitude cutoff. If the various redshift slices probe to different limiting masses 
(or luminosity) on the galaxy mass function, then the slice going to lower mass will have 
more objects per unit volume. This will obviously result in a larger pair fraction, i.e. a 
larger number of c ompanions within a specified radius, independent of any evolution (cf. 
Berrier et al.ll2006l ). For this reason we attempt to encompass the same "classes" of objects 
across the full redshift range of our study and do this by adopting a fixed cutoff in Ly. 
We point out in §5.3 that the use of a fixed cut in luminosity is more appropriate than 
a passive evolution cut (i.e. lower luminosity at low redshift). The vast majority of the 
merging systems are star-forming disk galaxies with ongoing star formation even to the 
lowest redshifts in our sample (§5.1); they are not passively evolving. 



3.2. Identification of GalcLxy Pairs 

To identify possible galaxy pairs, we searched for the nearest-neighbor galaxy for 
each galaxy in the sample - using photometric redshifts to limit the difference in redshift to 
Az < ±0.05 (~ the typical uncertainty in the photometric redshift estimates). Thus both 
the primary and its possible companion must have L > Ly and the most probable redshifts, 
from the photometric redshift fitting, must be within 5%. From computed probability distri- 
butions derived from the photometric redshift fitting, we estimate that ~90% of pairs would 
be found. This is likely to change slightly (ranging from ~80-95%) with redshift since the 
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error in the photometric redshift estimate depends on redshift. However, as discussed in §3.3, 
we calculate the expected random distribution using precisely the same numbers of galaxies 
in each redshift slice and therefore this small variation is corrected for. We also limited the 
range of projected separations to minimize contamination from chance, line-of-sight pairs. 

Galaxy pairs were identified in two ways - a primary search yieldi ng over 95% of 



the s ample was done using the ground-based photometric redshift catalog (IMobasher et al. 



20071): this was supplem ented with a search for very close pairs using the ACS-based catalog 



(ILeauthaud et al.l 120071 ). The ground-based catalog was the primary catalog used for our 
analysis since it contains the majority of the sample and immediately provides redshifts and 
absolute magnitudes for all objects. However, the resolution of the ground-based images 
mean that a small fraction of the most closely separated galaxy pairs are missed in the pho- 
tometric redshift-based search. To ensure that we have indentified all pairs at these close 
separations, we then used the catalog extracted from the ACS images to search for close 
companions which are blended and appear as a single entry in the photometric redshift cat- 
alog. (For these close objects, we adopted the single redshift from the photometric redshift 
catalog.) Any "ACS-only" pairs found in the secondary search which were already picked 
up in the ground based catalog were of course excluded. "ACS-only" pairs for which the 
photometric redshift fit was poor were also dropped from further consideration since these 
could be superpositions at different redshifts. Each of the "ACS-only" galaxy pairs was also 
visually inspected to remove "artificial pairs" which might have been present in the ACS 
catalog - for example, features of the same nearby galaxy, such as a bright spiral arm or 
residual cosmic rays. 



3.3. Correction for Line-of-Sight Superpositions 

Comparison of the observed pair distribution with that of a random sample is needed 
to understand whether the pair distribution could be simply due to chance line-of-sight 
coincidences. To evaluate the statistical significance of the measured distribution of nearest- 
neighbor galaxies, we generated samples of galaxies, randomly distributed over the field 
(excluding masked areas) but with redshift distribution identical to that of the observed 
galaxies. For ten repetitions of this Monte-Carlo simulation, we then searched for nearest 
neighbors at Az within ±0.05 out to a projected physical distance of 150 kpc, just as with 
the real galaxy sample. The number of galaxies in the Monte-Carlo simulations was adjusted 
to compensate for the variable volumes of the different redshift slices. For the Monte Carlo 
simulation, the spatial position of each real galaxy in the field was randomized but the 
redshifts were retained; thus any redshift dependent effects were automatically maintained 
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the same for the real galaxy sample and for the simulation galaxies. The ratio of nearest- 
neighbor distributions for the observed galaxies and the Monte-Carlo simulations is shown 
in Figure [T] for each redshift range. A completely random set of observed pairs would have 
a value of unity at all separations in Figure [TJ Instead, the distribution shows a clear 
excess of neighboring galaxies at small projected separations (and a deficiency at the larger 
separations). At close separations, the observed/random ratio is above unity and falls off to 
approximately match the random sample at separations larger than 40 — 50 kpc. 

It is important to recognize that since we are searching for the nearest (i.e. the first) 
neighbor to each L* galaxy, our measurements are sensitive only to close pairs and they 
are not sensitive to the large-scale structure as measured in the two-point correlation func- 
tion. We also note that the Monte-Carlo simulations were entirely random, i.e. without 
any superposed large scale clustering. Clearly some increases in the frequency of pairs are 
associated with the non-uniform large scale distribution or clustering of galaxies. But it 
would be wrong to artificially remove this component of the pair population if it is indeed 
affecting the actual pair fraction or merger rate. The large-scale structures probably do not 
have an important effect on the distribution of nearest neighbors as measured here, although 
obviously large structures will affect the large-scale correlation function. 

To minimize the contamination from "random pairs" (i.e. chance line-of-sight superpo- 
sitions), one should select only those pairs with apparent separations substantially less than 
40 kpc. For separations < 20 kpc, the expected contamination by random pairs is typically 
15% in each redshift bin (see Figured]). While we could decrease this percentage further by 
reducing the maximum separation, the total number of real pairs would then decrease, yield- 
ing a statistically less-significant sample. A minimum separation of 5 kpc was also adopted; 
at smaller separations, the ACS images cannot reliably discriminate cases where internal 
galaxy structures appear as multiple sources and hence an erroneous pair. 

4. RESULTS 

The results of our pair identification are summarized in Tabled! In total, 1,749 galaxy 
pairs were found across the entire COSMOS field. In Table dl we separate the numbers of 
pairs found in ground-based and ACS catalogs in each redshift bin and the expected random 
component, based on the Monte-Carlo simulations. The pair fraction as given in Tabled! is 
the fraction of galaxies with a nearest neighbor at 5 to 20 kpc, corrected for the expected 
random count. Since the pair fraction is calculated using the total number of galaxies in 
close pairs, in the case of a compact group, each galaxy is counted once, provided it is within 
20 kpc of another galaxy. Thus the pair fraction as given here accurately represents the 
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probability that an individual galaxy has a close companion. In cases where the nearest 
neighbor to one galaxy is a random projection, all galaxies within 20 kpc would still be 
picked up because they would each have a different nearest neighbor. This ensures that all 
pairings within 20 kpc are found and then we correct for chance projections as discussed in 
§3.3. 

In order to place our results in perspective, we have determined a low redshift galaxy 
pair fraction using the same criteria adopted for the COSMOS data. The COSMOS survey 
cannot be used for this "local" value since the volume sampled in the three lowest redshift 
bins (i.e. at z < 0.4) becomes too small while at the same time, evolution results in there 
being fewer pairs at low redshift. We A vere fortunate to b e able to use a published catalog of 
galaxy pairs produced from the SDSS (lAUam et al.ll2004l ). From the SDSS catalong, 45 pairs 
(out of a total of 705 cataloged pairs) at 2; < 0.1 met the selection criteria for the projected 
separations and luminosity (as were used for the COSMOS pairs). A spectroscopic redshift 
was available for at least one member of each of these pairs but not always for both members 
of the close pairs. Due to this low level of incompleteness, lAUam et al.l (120041 ) estimate that 
the expected contamination of the sample is < 3.4%. 

Figure [2] shows the galaxy pair fraction as a function of redshift for the COSMOS pairs 
together with the local value de t ermine d from the SDSS. The error bars shown represent 
la Poisson errors. IScoville et al.l (l2007al ) discuss the cosmic variance and Poisson noise for 
the COSMOS survey. It is seen there that the cosmic variance decreases relative to the 
Poisson noise at low masses and should not be significant on the scale of individual galactic 
halo masses. A simple least squares fit, weighted by the error bars, to the data in Figure [2] 
yields a slope of n = 3.1 ± 0.1. (If the SDSS point is excluded, the best fit power-law has 
n = 2.8 ± 0.1) We discuss this result below and compare our findings with previous work. 



5. DISCUSSION 

The results shown in Figure [2] imply strong evolution in the galaxy pair fraction of > Ly 
galaxies with redshift. The pair fraction rises by a factor of ~10.6 - from ~1% at z = 0.1 
to ~10% aX z = 1.1 - and if extrapolated to higher redshift, would imply that nearly half of 
all > Ly galaxies at 2; ~ 2 are in close pairs (!) 

The power-law fit shown in Figure [2] provides a reasonable representation of the data, 
with the majority of the data points lying within 2a of the straight line of slope n = 3.1. 
The only really significant deviation is for the data point at z = 0.35; this turns out to 
be the redshift bin corresponding to several major galaxy clusters in the COSMOS field 
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( IScoville et al.ll2007af ). Thus this high data point is probably high due to an increased pair 
fraction associated with the overdensity of galaxies in these extended clusters. Although we 
discuss the spatial distribution of the identified pairs below, a more thorough investigation 
of the pair distribution with respect to local environment needs to be carried out to properly 
address this question; this will be the subject of a second paper. 

It is important to note that had our study been limited to a smaller range of redshift, 
smaller area, or shallower depth, the value of n might have been significantly different due 
to cosmic variance and Poisson statistics. For example, our sample, limited to 2; < 0.5, 
would have found n ~ 5 while the intermediate redshifts z = 0.3 — 0.8 would have yielded 
n ^ 0. Reduced volumes and more limited redshift ranges can be severely affected by 
counting statistics and large-scale structure excesses or deficiencies (i.e. cosmic variance). 
The COSMOS survey provides superb angular resolution and depth with the COSMOS ACS 
and deep ground-based images with excellent seeing over the very large field. In addition, 
the availability of a robust SDSS data point allows us to pin down the local value for the 
pair fraction. 



5.1. Morphology and Spectroscopy of Pairs 

To validate the sample of galaxy pairs identified here, we visually inspected the ACS 
images of all pairs. Figure H] shows the ACS images for a sample of 25 galaxy pairs at a variety 
of separations ordered by their redshifts (but selected at random). The projected separation 
along with the mean photometric redshift for the pair is shown in the lower right corner of 
each panel. Similar plots were made in order to visually inspect all of the 1,749 galaxy pairs in 
the COSMOS field. These plots were initially made for the purpose of screening out spurious 
detections (due to residual cosmic rays, etc.); however, the depth and resolution of the ACS 
images also makes them extremely useful for discussing the morphology of the galaxies. A 
substantial fraction of these pairs are obviously interacting, showing clear signs of disturbed 
morphologies (tidal bridges, tails, etc.). Their morphologies are indicative of early/mid- 
stage major mergers, involving relatively equal mass pairs, as judged by comparison with 
the morphologies of strongly interacting pairs at low redshifts. 

We are also in the process of obtaining spectra for the complete sample of pairs. Optical 
spectra for a small subsample of pairs have recently been obtained using DEIMOS on Keckll. 
A preliminary analysis of these data indicates that >85% of our photometrically identified 
pairs are indeed galaxies with AV < 600kms~^, which is in agreement with our expectation 
based on the expected numbers of observed vs. random pairs shown in Figure [H (The 
complete catalog of objects and their images and spectra will be made available electronically 
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in a future paper.) 



5.2. Spatial Distribution of Pairs 

Figure [3] shows the distribution of all galaxy pairs found in the COSMOS field in each 
of the /S.Z = 0.1 redshift bins. The distributions of pairs are clearly not random across the 
field. The most massive large-scale structure seen in the COSMOS field at z ~ 0.73 lies in 
the 0.7 < z < 0.8 redshift bin and a higher density of pairs can be seen in the northwest 
of the field in this panel in Figure [3l Similarly, there is a significant excess of pairs in the 
z = 1.05 slice on the west side of the field. At the higher redshifts {z > 0.5), the COSMOS 
field contains very large volumes in each redshift slice so that the overall contribution of a 
single structure or cluster is diluted and contributes < 10% to the total counts; at lower 
redshifts, an increase in the pair fraction within large clusters (or alternatively a decrease 
in the pair fraction in large voids) can cause a larger deviation in Figure [2] simply due to 
the smaller comoving volume of the field. This could partially account for the higher count 
seen in the z ~ 0.35 bin (Figure [2]). However, the overall number of pairs in the low z 
slices is also low (due to both the small comoving volume and reduced pair fraction at low 
redshift), hence, the deviations can also be largely due to small number statistics. In order 
to assess the significance of cosmic variance across the field, we attempted to divide the field 
into separate quadrants which were then individually analyzed. Unfortunately, the results in 
the separate quadrants were even more strongly affected by small number statistics and the 
quadrant analysis is even less meaningful in assessing cosmic variance. For now, we simply 
note that one may need to sample volumes at least as as large as 10^ Mpc^ in order to average 
out the variations due to LSS, something that COSMOS was designed to do at ^ > 0.5, as 
the SDSS does at 2; < 0.1. Nearly the same slope (n ~ 3.1) would have been found for the 
data in Figure [2] had only the SDSS data point and the COSMOS data at z > 0.5 been used 
in the fit. 



5.3. Comparison with Previous Studies 



To facilitate comparison of our results for evolution of the galaxy pair fraction (oc 
(1 + ^)3-i±o i) with prior studies, we compile in Table [2] a summary of earlier studies of 
pair fraction and merger rate evolution. This table also includes several recent studies which 
used morphological criteria presumed to be related to interactions and merge rs as a basis for 



direc t ly sampling systems that have interacted and subsequently merged (e.g.. lLe Fevre et al. 



2OOOI : IConsehce et al.ll2003l : iLaverv et al.ll2004J : iLotz et al.ll2006f ). A critique of both methods 
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(pair fraction and disturbed morphologies) is provided in lLavery et al.l (120041 ) . The results of 
both are usually parameterized as a power-law in (l + z). Although it is reasonable to believe 
that the two power laws would be related, the indices are probably not the same since not 
all pairs will merge and their merging timescales are different. The fraction of pairs which 
do eventually merge must depend on their separations and the distribution of separations 
is likely to vary with redshift due to merging. Several studies have argued that power-law 
index m for the m erger rate ranges be tween ±1 from the index n used to parameterize the 



pair fraction (e.g., iLavery et al.ll2004l ). For samples like ours, involving pairs of luminous 
(i.e. > Ly) and presumably massive galaxies at dose separations (< 20kpc), we expect that 
(m — ra) ~ since most of these systems should eventually merge. 

Most of the investigations listed in Tables [2] used spectroscopic redshifts (at least for 
objects at 2; < 1) so the samples are relatively small - typically ~ 2 — 8% the size of our 
COSMOS sample. These samples may also be incomplete at high redshift where spectro- 
scopic surveys are unable to resolve very close pairs. While some of these studies suffer 
serious incompleteness due to their apparent magnitude limit, several of the most recent, 
using deep surveys, have sufficient sensitivity to adopt an absolute magnitude selection limit 
similar to what we have done here. The derived values of m and n span an enormous range 
- from — 6. While "m versus n" can account for some (i.e. ±1) of this variation, it is 
more likely that different redshift ranges, different ways of correcting for non-physical (line of 
sight) galaxy pairs, cosmic variance, and small sample sizes are the major contributors to the 
large discrepancies in power-law indices. The straight average of all measurements in Table 
12] (excluding the current work) is ~ 2.1 ±0.8, with the mean for methods using pair-fraction 
and for methods using morphology being nearly equal. One notable difference concerns 
those studies that have adopted "passive luminosity evolution" (PLE) models to "correct" 
the observable luminosity measurements to mass measurements. These more recent studies 
(flagged in Table [2]) give systematically smaller values, with an average of 1.0 ± 0.3. 



5.4. Galcixy Masses and Luminosities: Passive Luminosity Evolution (PLE)? 



Many recent studies of the merger rate have assumed passive (or pure) luminosity evo- 
lution (PLE) in an attempt to select similar galaxies at all redshifts. PLE is typically 
parameterized by r nagnitudes M(z) = M(z = 0) — Qz where the empirical assumption has 
been that Q — 1 ( Lin et al. 2004 : Patton et al. 2OO2I ). We tested such a model, selecting 
galaxies using My = —19.8 in our highest redshift bin (1.1 < z < 1.2), but including fainter 
galaxies at the lower redshifts, assuming Q = 1 in the equation above. The rest of the 
selection criteria remained the same, resulting in a sample of 117,015 galaxies. (The sightly 
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increased number counts arise in the lowest redshift bins.) Using this PLE sample of galaxies 
and counting the pairs again, we obtained a fit with n = 2.2 ± 0.1. 

However, we do not believe that PLE is justified empirically. To test for evolutionary 
variation in the ratio M/Ly as a function of redshift, we made use of the galaxy masses 
derived from the CO SMOS photometry. To e stimate the mass of each galaxy in the COSMOS 
photometric catalog, iMobasher et al.l f|2007n used a relation between M/ Ly and rest-frame 
{B — V)q colors given by iBell et al.l (120051 ). These masses are shown for the pair sample 
galaxies in Table [3l They are consistent with relatively little or no variation in mean mass 
versus redshift for our pair sample. The median spectral type of the galaxy pairs indicates 
that they are late-type spirals for which passive evolution is inappropriate since these late 
type galaxies have continuous, on-going star formation. The apparent difference in slope is 
more likely due to an increased merger fraction at faint magnitudes than passive evolution for 
the massive galaxy pairs. In summary, we find no basis for adopting PLE for the COSMOS 
sample of galaxy pairs and instead, retain the fixed Ly sample as most relevant for analysis 
of the evolution of galaxy merging. 



6. SUMMARY 

COSMOS HST and ground-based images covering the full 2 0° COSMOS field have 
been used to identify an unprecedented sample of 1,749 luminous (> Ly) galaxy pairs at 
redshift z < 1.2. These data, supplemented with a complete sample of 45 "local" galaxy 
pairs selected from the SDSS, have been used to determine the pair fraction of luminous 
galaxy pairs as a function of redshift over the range < z < 1.2. In addition, galaxy masses 
obtained from an analysis of the extensive multi-wavelenghth (UV-opt-NIR) COSMOS data, 
have been used to investigate possible evolution in M/L versus redshift. From our analysis, 
the following conclusions are drawn: 

1. Evolution of the galaxy pair fraction of close pairs (with projected separation of 5 — 20 
kpc) can be fit as (1 + 2;)", where we find n = 3.1 ± 0.1. 

2. The HST-ACS images suggest that a large majority of our luminous galaxy pairs are 
strongly interacting/merging disks galaxies. 

3. The mean M/L ratio of the galaxy pairs in our sample shows no evidence of evolution 
with redshift. This argues against the use of passive evolution models in conjunction 
with the sample selection. (In any case, it is hard to justify using passive evolution 
models for star forming disk galaxies which are typically the focus of merger galaxy 
studies). 
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The high angular resolution and large area coverage provided by the COSMOS HST-ACS 
imaging were crucial in allowing us to cleanly separate close pairs (separations < 20 kpc) out 
to z = 1.2 and providing an enormous sample of 1749 galaxy pairs - 10 — 20 x larger that that 
used in previous studies. The ability to focus the study on close pairs minimized the number 
of chance line-of-sight pairs and increased our ability to identify strongly interacting/ merging 
systems. The very large sample size was critical for reducing cosmic variance and Poisson 
errors; these variances have been a major problem for earlier studies. 

The resolution of the ACS images would, in principle, allow us to extend the investiga- 
tion to 2; ~ 2 once deeper near-infrared (F, J, H, Kg) and Spitzer-IRAC (3.6, 4.5, 5.6, 8.0 fim) 
data are available for the photometric redshift analysis. The results presented here show that 
~ 10% of all luminous (> Ly) galaxies are in close pairs already by 2; = 1.2; extrapolating 
the (1 -I- z)^-^ power-law to higher redshift suggests that ~50% of all luminous galaxies are 
in close pairs by 2; = 2. Thus major interactions and mergers must represent an important 
and very likely dominant process of galaxy evolution in the early universe. 

The HST COSMOS Treasury program was supported through NASA grant HST-GO- 
09822. We gratefully acknowledge the contributions of the entire COSMOS collaboration 
consisting of more than 70 scientists. We also wish to thank the referee for helpful comments 
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thank Lisa Kcwlcy for many helpful discussions concerning this project as well as Sahar 
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would like to thank the Institute for Astronomy, University of Hawaii for support during a 
sabbatical stay where this project began. DBS and NZS also acknowledge the hospitality of 
the Aspen Center for Physics. 
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Table 1. Number of COSMOS Galaxies (> Ly) with a Companion at Projected 

Separation 5 — 20 kpc 





0.15 


0.25 


0.35 


0.45 


0.55 


0.65 


0.75 


0.85 


0.95 


1.05 


1.15 


ACS only 

















5 


9 


14 


44 


85 


59 


Ground Based 


10 


27 


148 


117 


163 


381 


597 


398 


757 


772 


837 


Total 


10 


27 


148 


117 


163 


386 


606 


412 


801 


857 


896 


Total - Random 


8 


19 


117 


93 


124 


244 


424 


299 


584 


612 


642 


Galaxy Sample 


386 


1043 


2956 


3097 


3680 


7320 


9575 


6719 


8888 


9276 


6281 


Pair Fraction 


0.021 


0.018 


0.040 


0.030 


0.034 


0.033 


0.044 


0.045 


0.066 


0.066 


0.102 


Gomoving Volume sampled (10® Mpc^) 


0.09 


0.23 


0.41 


0.61 


0.81 


1.01 


1.21 


1.38 


1.54 


1.69 


1.82 



Gentroid of redshift bin of width 0.1 
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Table 2. Summary and Comparison with Previous Work 



Reference mag^ z Range'^ Sample Size'^ #P or M° m or n 



Galaxy Pair Fraction Studies 



Zepf & Koo (1989) 


B <22 


A 


< z>= 0.25 


~ 1000 


20 


4.0 ±2.5 


Carlberg et al. (1994) 


V < 22.5 


S 


< z >= 0.42 


1062 


14 


2.3 ± 1.0 


Burkey et al. (1994) 


7=18.5-23 


A,S 


0.4-0.7 


146 


50 


3.5 ±0.5 


Woods et al. (1995) 


/ < 24 


N 




~ 1000 


23 





Neuschaefer et al. (1995) 


/ < 22 




^med — 0-5 


~ 4500 







Yee & EUingson (1995) 


r < 21.5 


S 


<z>= 0.38 


107 


25 


4.0 ± 1.5 


Neuschaefer et al. (1997) 


7=18-23 


S 


<l-2 


~ 22400 


90 


1.2 ±0.4 


Patton et al. (1997) 


r < 22 


s 


< z>= 0.33 


545 


73 


2.8 ±0.9 


Wu & Keel (1998) 


/ < 24 


A 


~ 2 




10 


0-2 


Carlberg et al. (2000) 


Mr < 19.8 


s 


0.1 - 1.1 


~3300 


109 


0.1 ±0.5 s 


Le Fevre et al. (2000) 


5m < 1.5 


s 


< 1 


285 


26 


2.7 ±0.6 


Patton et al. (2002) 


Rc < 21.5 


s 


0.12 - 0.55 


4184 


88 


2.3 ± 0.7s 


Lin et al. (2004) 


Ms=-(21-19) 


s 


0.45 - 1.2 


5000 




1.60 ±0.29 


Lin et al. (2004) 












0.51 ±0.28 


This work 


> Lv 


Ph 


1.2 


106118 


1749 


3.1 ±0.1 



Morphologically Selected Merger Studies 



Le Fevre et al. (2000) 
Reshetnikov (2000) 
Conselice et al. (2003) 
Lavery et al. (2004) 
Lotz et al. (2006) 



5m < 1.5 



Mb < -18 



> 0.4Lt 



< 1 
0.5 - 1.0 
0-3 
Ph,S 0.1-1.0 
S 0.2-1.2 



S 

Ph 
Ph,S 



285 
1212 
2368 



49 
14 
43 
25 
157 



3.4 ±0.6 

3 6+^-2 
■^■"-0.9 

4-6 
5.2 ±0.7 
1.12±0.60S 



^Magnitude selection limits as stated in the reference. 5m refers to difference in magnitude and Rc is 
a Kron-Cousins R magnitude. 

''Method used for redshift measurement: S = spectroscopic, Ph = photometric, A = assumed, N = 
none 

'^Redshift range as stated in the reference. In some cases, only the mean, median, or upper limit of 
the redshift range was given and then compared to a local value. 

"^Galaxy sample size, as stated in the reference. In the case of lLaverv et al. I l l2004h . 162 WFPC2 fields 
were used. 

"^Total number of identified Pairs (P) or morphologically selected Mergers (M) in the sample. 

'Exponent: merger rate (m) or pair fraction (n), assuming evolution oc (1 + z)™'". Error bars are 
assumed to be Icr. 

^Assuming passive luminosity evolution. 

''Values for total number of galaxies and pairs includes only the COSMOS data. Including the SDSS 
data used for the z = — 0.1 data point would only increase these values 



Table 3. Log of Stellar Masses of Galaxy Pairs 



z^- 0.15 0.25 0.35 0.45 0.55 0.65 0.75 0.85 0.95 1.05 1.15 



Mean 10.03 10.04 10.07 10.11 10.15 10.14 10.11 10.15 10.13 10.13 10.09 
Median 9.83 9.98 10.02 10.11 10.17 10.11 10.09 10.13 10.12 10.12 10.07 
Stddev 0.48 0.43 0.45 0.42 0.47 0.45 0.46 0.44 0.44 0.46 0.45 



''Centroid of redshift bin of width 0.1 
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Fig. 1. — Ratio of observed to random pairs versus pair projected separation in Az — 
0.1 redshift bins between z = 0.1 and 1.2 measured for galaxies from the ground-based 
photometric redshift catalog. This distribution represents the likelihood that a galaxy has 
its nearest neighbor at a given separation referenced to the likelihood obtained for the same 
number of galaxies randomly distributed over the same area. Since we are searching for the 
nearest (i.e. the first) neighbor to each L* galaxy, these measures are sensitive only to close 
pairs, not to the large-scale structure as measured in the two-point correlation function. 
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Fig. 2. — Fraction of > Ly galaxies in close pairs (5 — 20 kpc) versus (l + z) for the COSMOS 
field. The star marks the local {z = — 0.1) data point determined using data from the 
SDSS (see text). Vertical bars represent la error. The straight line least squares fit to the 
data has a slope of n = 3.1 ± 0.1. 
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Fig. 3. — Distribution of galaxy pairs across the entire 2 0° COSMOS field in Az — 0.1 
redshift bins. The + symbols represents the centroid position of each identified pair. 
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Fig. 4. — ACS I-band (F814W) postage stamps (10" x 10") of a representative sample of 
25 luminous (> Ly) galaxy pairs with projected separations < 20 kpc in the COSMOS 
field ordered by their mean photometric redshift. The redshift of the pair and projected 
separations in kpc are shown in the lower right corner of each panel. The ACS source 
positions for each member of the pair are marked by green crosses. 



